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ABSTRACT 

We report the first space-based very deep near-infrared continuum observations of a uniform sample 
of z = 6.5 galaxies with log(L(La)) = 42.5 — 43 erg s _1 selected from narrow-band line searches and 
with spectroscopically confirmed Lya emission. The 1.4/im HST WFC3 observations are deep enough 
(AB(ler) = 28.75) to measure individual continuum magnitudes at this redshift for all of the objects. 
We compare the results with continuum-selected samples at the same redshift and find that Lya 
emission is present in 24% of all galaxies with Mab(L350A) < —20 at z = 6.5. The error in this 
quantity is dominated by systematic uncertainties, which could be as large as multiplicative factors 
of three. The Lya galaxies are extended but small (size < 1 kpc), and have star formation rates of 
approximately 10 Mq yr~ x . We find a mean L(hya)/vL u at 1400 A of 0.08, with the seven objects 
showing a range from 0.026 to 0.26, implying that there is little sign of destruction of the Lya line. 
All of the properties of the z = 6.5 sample appear to be very similar to those of Lya emitters at lower 
redshifts. 

Subject headings: cosmology: observations — early universe — galaxies: evolution — galaxies: for- 
mation — galaxies: high redshift 



1. INTRODUCTION 

At all redshifts from z = 0toz = 7a fraction of 
galaxies can be detected in the Lyman alpha line of hy- 
drogen (e.g. Cowie and Hu 1998; Deharveng et al. 2008; 
Hu et al. 2010; Ouchi et al. 2010; Vanzella et al. 2010). 
The Lya line is easily detected in narrow-band imag- 
ing or blind spectroscopic surveys, and Lya emission- 
line searches have been widely used to find high redshift 
galaxies. In addition, for the highest redshift galaxies, 
this line is the only spectroscopic signature that can be 
used to confirm the redshift of a galaxy selected on the 
basis of color properties. However, Lya is a difficult line 
to interpret. Because the line is resonantly scattered by 
neutral hydrogen, determining its escape path and there- 
fore its dust destruction is an extremely complex prob- 
lem both theoretically (e.g., Neufcld 1991; Finkelstein et 
al. 2007) and observationally (e.g., Kunth et al. 2003; 
Schaerer & Verhamme 2008, Ostlin et al. 2009). It is 
therefore critical to determine the ultraviolet continuum 
properties of the Lya emitting galaxies (LAEs) in order 
to understand how to interpret them in terms of physical 
quantities such as the star-formation rate and to relate 
them to UV-continuum selected samples at the same red- 
shift (e.g., Bouwens et al. 2010). 

Extending this measurement to the highest redshifts 
possible is of considerable interest because the LAEs may 
be used as probes of the structure of the intergalactic 
gas at z > 6 (e.g., Dijkstra & Wyithe 2010, Tilvi et 
al. 2010, Dayal et al. 2011, Laursen et al. 2011). Since 
Lya photons resonantly scatter off of neutral hydrogen, 
the Lya emission from galaxies scatters from the line of 
sight as it propagates through the intergalactic medium 
(e.g. Loeb & Rybicki 1999; Zheng et al. 2010) potentially 
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reducing the strength of the observed Lya. Thus the 
frequency of Lya emitters relative to the UV continuum- 
selected galaxy population is potentially a function of 
the degree of ionization of the intergalactic medium and 
could be used to probe the epoch of reionization of the 
universe, at redshift z > 6 (e.g. Stark et al. 2010, 2011). 

Here we extend the measurement of the UV-continuum 
properties of the LAEs to z = 6.5, the highest redshift 
where substantial samples of these objects exist, using 
very deep near-infrared continuum observations of a uni- 
form sample of z = 6.5 LAEs obtained with the WFC3 
camera on HST. These are the first observations deep 
enough to measure individual continuum magnitudes at 
this redshift for all the objects in this population and we 
find that their UV-continuum properties and the frac- 
tion of Lya emitters relative to UV-continuum selected 
samples are very similar to that seen at lower redshift. 
The observations show no indication that the onset of 
reionization is at z = 6.5. 

2. OBSERVATIONS 

We obtained 1.4 /im (rest frame 1850A) continuum 
detections for a sample of 7 Lya emitters (LAEs) at 
z = 6.5 in a 600 arcmin 2 region centered on the mas- 
sive z = 0.3 galaxy cluster Abell 370. The sample was 
drawn from Hu et al.'s (2010) Lya spectroscopic atlas of 
88 z — 5.7 and 30 z = 6.5 galaxies. The Hu et al. Lya 
emitters were chosen based on their narrow-band excess 
and their broad-band colors in very deep images obtained 
with the Suprime camera on the Subaru 8m telescope. 
They were then spectroscopically confirmed using the 
DEIMOS spectrograph on the Keck 10m telescope. They 
form a sample of objects with Lya rest frame equivalent 
widths above 20A and log(L(La))= 42.5 - 43 erg s 1 in 
the selected redshift interval. These objects have narrow 
Lya emission lines (velocity widths 150 - 300 km s _1 ) 
and, as we shall show below, while they are small, they 
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are marginally spatially extended in the present obser- 
vations. We therefore assume that they are powered by 
star formation and not by an active galactic nucleus. 

We used the WFC3 camera on HST with the F140W 
filter to obtain 1.4^m observations of all of the z = 6.5 
objects in the A370 held of that survey. Each of the 
z = 6.5 galaxies was observed for 2 orbits using the 
F140W filter of the WFC3 camera on HST (Mackenty 
et al. 2010). A four point dither pattern was used 
and repeated 3 times, giving a total exposure time of 
4235s for each field. We used the combined drizzled im- 
ages together with the calibration produced by the HST 
pipeline, but then removed a further background deter- 
mined from the median average of the fields. The mea- 
sured FWHM for faint stars in the images is 0.22 arcsec- 
onds. As we shall discuss below, the galaxies are small 
so, given the high spatial resolution of the HST images, 
we used 0.8" diameter magnitudes and corrected these 
to approximate total magnitudes using an offset of 0.15 
magnitudes computed from the enclosed energy curves 
of the instrument. We measured the noise in the images 
by measuring the fluxes in random blank sky positions 
and found ale noise level of 28.75 for the corrected 0.8" 
diameter magnitudes. 




Fig. 1. — Comparison of the HST F140W observation for atypical 
object (HC023939-013451) with the ground based data from Hu et 
al. (2010). The lower left panel shows the F140W image, the lower 
right the 9120A narrow band image, the upper left the z band 
image and the upper right the R band image. The circle marking 
the LAE is one arcsecond in radius. 



All of the galaxies in the sample were detected; a typ- 
ical example is shown in Figure [TJ The faintness of the 
image, even in a very deep z band exposure such as that 
shown (5cr ~ 25.4 AB), illustrates the difficulty of mak- 
ing accurate continuum measuruments from the ground. 
The observed continuum magnitudes and other quanti- 
ties are listed in Table [TJ The first detected z — 6.5 
galaxy, HCM6A (Hu et al. 2002), is a member of the 
sample and since this galaxy lies in the strong-lensing 
region of A370, for this object we also give the values 
corrected for lensing amplification (quantities in paren- 



theses in Table [JJ . All of the remaining galaxies lie far 
outside the strong lensing region and have negligible am- 
plification. 

For a comparison sample at z = 5.7 chosen in the same 
way and with comparably deep continuum magnitudes, 
we took Hu et al.'s sample of six z = 5.7 spectroscopically 
confirmed objects in the GOODS-N field, using published 
F850LP data obtained with the ACS on the HST (Gi- 
avalisco et al. 2004) to define our continuum magnitudes. 
The z = 5.7 magnitudes and errors are the auto mags 
from the catalog of Giavalisco et al. (2004); they also pro- 
vide a good representation of the total magnitude. The 
rest-frame wavelength corresponding to the pivot wave- 
length of this filter is 1350A. While the F140W filter has 
no response at the Lya line, the F850LP filter has 10% 
of peak response at 8140 A (Lya for z — 5.7). However, 
the emission-line contribution to the continuum fluxes 
is small (< 20%) even for the highest equivalent width 
(EW 140A) z — 5.7 galaxies. The properties of these 
objects are also listed in Table [TJ 

3. DISCUSSION 

The LAE galaxies at both z — 5.7 and z = 6.5 are 
extremely compact (Figure [5]) . For each redshift sample 
we formed a stack by summing all the images, center- 
ing on the peak flux. For the z = 6.5 objects we ex- 
cluded HCM6A, which is lensed, and HG023949-013121, 
which lies close to a bright galaxy. A comparison of the 
stacked continuum images of the two samples with the 
point spread function (PSF) determined from faint spec- 
troscopically confirmed stars in the two fields, shows that 
images at both redshifts are extended: the deconvolved 
FWHM is 0.08" at z = 5.7 and 0.14" at z = 6.5, cor- 
responding to sizes between 0.5 and 1 kpc. This is very 
similar to values measured for Lya emitters (LAEs) at 
z = 3.1 (Bond et al. 2009, 2010), where both the UV 
continuum and the Lya emission are seen to be sub-kpe. 
The sizes are smaller than those found by Taniguchi et al. 
(2009) at z = 5.7 but the Taniguchi measurements were 
made in a bandpass containing the Lya line and it is 
possible that this could be increasing the measured size 
in their sample if the Lya emission is more extended. 
The UV size measurements provide a reference for fu- 
ture measurements of the Lya sizes, which will allow 
us to measure the extent of the Lya scattering region 
in galaxies. Such measurements are difficult but may 
be possible with ground-based adaptive optics or space- 
based narrow-band measurements. 

Figure [3] shows our primary result. For both samples 
we plot the ratio of Lya luminosity to continuum lu- 
minosity redward of Lya (or equivalently the rest-frame 
equivalent width of the Lya line, to which the ratio is ap- 
proximately proportional [Dijkstra & Westra 2010]) as a 
function of Lya luminosity. Following Stark et al. (2010), 
we have assumed f v ~ v 3 in computing the equivalent 
width. Adopting f v ~ i/~ - 8 as suggested in Taniguchi et 
al. (2009) would increase the equivalent widths by a mul- 
tiplicative factor of 1.4 for the z = 6.5 systems and 1.1 for 
the z = 5.7 systems. However, we note that we focus less 
on the Lya equivalent width than on the UV continuum 
luminosity, which is less sensitive to the assumed /„ since 
the wavelength difference is smaller. Here too we assume 
/„ ~ v a in comparing the two redshifts. The z = 5.7 
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TABLE 1 
Galaxies in the sample 



Name 


RA(2000) 


Dec(2000) 


Redshift 


logL(Lya) 


mag a 


Rest EW 6 










(erg s" 1 ) 


(AB) 


(A) 


HC023927-013523 


39.863293 


-1.589833 


6.450 


42.73 


26.68 ±0.16 


64 ±9 


HC023939-013432 


39.914417 


-1.575667 


6.549 


42.76 


26.03 ± 0.09 


36 ±3 


HC023939-013451 


39.914791 


-1.581028 


6.531 


42.90 


26.73 ±0.17 


97 ± 15 


HC023949-013121 


39.956711 


-1.522667 


6.564 


42.99 


28.00±°;S 


otr 7 +329 

•"'-lie 


HCM6A 


39.978001 


-1.559111 


6.559 


43.41(42.75 c ) 


24.62(26. 27 c ) ±0.02 


44 ±1 


HC024001-014100 


40.007500 


-1.683388 


6.545 


42.72 


27.84±°;« 


175±« 4 


HC024004-012252 


40.017708 


-1.381278 


6.502 


42.62 


26.37 ±0.12 


37 ±4 



HC123607±620838 
HC123558+621017 
HC123613±620748 
HC123651±621936 
HC123652+622152 
HC123717±621759 



189.033004 
189.045471 
189.056106 
189.215240 
189.216751 
189.324677 



62.14394 
62.17144 
62.12994 
62.32683 
62.36460 
62.29974 



5.640 
5.672 
5.635 
5.675 
5.689 
5.663 



42.95 
42.58 
42.95 
42.96 
42.92 
42.91 



26.36 ±0.12 
26.54 ±0.18 
25.99 ±0.15 
26.89 ±0.24 

26.37 ±0.22 
26.52 ±0.15 



86 ±10 
42 ±8 
60 ±9 
140 ± 34 
79 ±18 
89 ±13 



Note. — (a) Continuum magnitude in the AB system at 1.4 fim (z = 6.5) and 0.91 [an (z = 5.7). (b) Rest-frame equivalent width of 
the Lya line calculated assuming a fiat-f,, spectrum for the continuum at wavelengths between Lya and ~ 1850A (rest-frame; z = 6.5) or 
~ 1350A (rest-frame; z = 5.7). (c) Values corrected for lensing amplification. 
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Fig. 2. — Stacked rest-frame near-UV images for (panel a) the 
z = 5.7 Lya-selected galaxies and (panel b) the subset of isolated 
and unlensed z = 6.5 Lya-selected galaxies. The point spread 
function measured from faint stars in the field is shown in the 
lower right corner. In both cases the images are resolved, though 
the effect is much more easily seen in the higher resolution ACS 
images of the z = 5.7 galaxies. 
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Fig. 3. — Fraction of light coming out in the Lya line versus 
that coming out in UV continuum redward of Lya. (The latter is 
calculated assuming a flat-/,/ continuum). This fraction is linearly 
related to the Lya equivalent width which is given on the right 
hand vertical axis of the figure. Objects at z = 6.5 are shown with 
blue diamonds and objects at z = 5.7 with red squares. The error 
bars are ±1 a. 



sample has identical median and mean values of 0.07 for 
the Lya to UV continuum luminosity ratio, while the 
z = 6.5 sample has a median value of 0.05 and a mean 
value of 0.08. Both are very similar to that expected 
from a galaxy undergoing constant star formation. For a 
Salpeter initial mass function, the latter value is ~ 0.065 
almost irrespective of metallicity (Chariot & Fall 1993; 
Schaerer & Verhamme 2008). This value is shown as the 
solid black line in Figure [3l There is little sign of evolu- 
tion with redshift, therefore, and, for this sample, little 
sign of destruction of the line relative to the continuum. 
The observed values catter above and below the aver- 
age values by multiplicative values of roughly three. The 
high values may correspond to younger age, top-heavy 
IMF or particular escape geometry (e.g. Finkelstein et 
al. 2007). 



The high value of the Lya to UV continuum ratio is 
likely to be, at least in part, a selection effect: if there 
is a range in the amount of destruction of Lya relative 
to the continuum, then the most luminous Lya emitters 
will be those with the least destruction, which may then 
have near-standard ratios. In particular, objects with- 
out Lya detections at the rest-frame EW limit of 20 A 
should not be analyzed with this result, which holds only 
for objects in the LAE samples. However, the result does 
imply that many objects with little Lya destruction do 
exist and that, for these luminous Lya-selected objects, 
we can convert between Lya and UV continuum lumi- 
nosities. At hrst sight this may seem inconsistent with 
the blue asymmetric profiles seen in these Lya emitters 
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(Hu et al. 2010). However, these profiles may be intrin- 
sic to the galaxy and produced by scattering in outflows 
(Verhamme et al. 2006) rather than being a consequence 
of Lya scattering in the surrounding intergalactic gas. 

Our measured ratio of Lya luminosity to UV con- 
tinuum luminosity can be combined with the 'con- 
ventional' Kennicutt (1998) relation between the star 
formation rate and the UV luminosity to obtain 
log SFR(M Q yr _1 ) = logL(Lya) - 42.03, where L(Lya) 
is in erg s _1 and, again, the value is for a Salpctcr initial 
mass function. The effects of metallicity may change this 
relation by as much as 0.2 dex (Schaerer & Verhamme 
2008) and it is also dependent on the exact star formation 
history. However, it can be seen from the Lya luminosi- 
ties in Tabled] that typical star formation rates, with no 
extinction correction, are in the range of 5-10 Mq yr _1 , 
where we have used the relationship above in computing 
the SFRs. The UV continuum-selected galaxies with- 
out detected Lya may have much more substantial Lya 
destruction and for these the conversion rate from Lya 
luminosity to SFR may be larger than the relation given 
above. 




-23 -22 -21 -20 -19 -18 



M w (1350A) 

Fig. 4. — Lya-selected 1350A continuum luminosity functions at 
z = 5.7 (red squares) and z = 6.5 (blue diamonds) derived from the 
Hu et al. (2010) Lyc* sample at these redshifts. Errors are ±ltr and 
upper limits are lc and shown with a downward pointing arrow. 
The black filled triangles show the z = 6.8 luminosity function de- 
rived in Bouwens et al. (2010) and the open downward pointing 
triangles that of Ouchi et al. (2009) while the solid black curve 
shows the Bouwens et al. (2010) maximum likelihood Schechtcr- 
function fit to the combination of these two samples. The blue 
dashed curve is this fit renormalized by a factor of 0.24 in number 
density to match the Lya-selected z = 6.5 continuum luminosity 
function. The red dotted curve is the z = 5.9 color-selected contin- 
uum luminosity function from Bouwens et al. (2010) renormalized 
by a factor of 0.18 in number density to match the z = 5.7 Lya- 
sclected continuum luminosity function. 



We can also now test whether the Lya emitters are 
comparable in continuum luminosity to samples obtained 
using continuum color-selection. From the full LAE sam- 
ple of Hu et al. (2010) we constructed continuum lu- 
minosity functions of the Lya-selected sample at both 
z — 5.7 and z = 6.5, using the accessible volumes com- 
puted for each object in Hu et al. We converted Lya 
luminosity to a UV continuum luminosity at 1350A (rest- 



frame) using the constant star-formation value of 0.065 
for the ratio of Lya luminosity to UV continuum light 
redward of Lya, and assuming the continuua are roughly 
constant in f v from 1350A to 1850A. These luminosity 
functions are shown in Figure 2] where we compare them 
with continuum luminosity functions determined from 
color-selected samples (Ouchi et al. 2009; Bouwens et al. 
2010). The Lya-selected objects have high continuum 
luminosities (or equivalently star formation rates), com- 
parable to the brightest continuum-selected objects at 
these redshifts. Our UV continuum luminosity function 
is lower by a factor of 5 than those derived by Shimasaku 
et al. (2006) and Kashikawa et al. (2006) whose results 
would imply that all UV continuum-selected objects at 
this redshift would also be LAEs. Both of these results 
may be due in part to the extreme difficulty of measur- 
ing accurate continuum magnitudes from ground-based 
data, particularly when the Lya lines lie within the con- 
tinuum bandpass. However, the difference is also due in 
part to the high Lya luminosity functions calculated in 
these papers, which are about a factor of two above that 
of Hu et al. (2010). Hu et al. (2010) argue that the differ- 
ence between their Lya LF and the previous determina- 
tions is a consequence of using a fully spectroscopically 
confirmed sample which eliminates some of the photo- 
metrically selected objects but Kashikawa et al. (2011) 
argue against this interpretation. It is unclear therefore 
which is the better Lya LF to choose, but we note that 
using the higher Kashikawa et al. (2011) Lya LF would 
result in a very large fraction of UV continuum-selected 
galaxies having Lya emission at near the case B ratio. 
This would be a rather surprising result, given that we 
expect at least some destruction owing to scattering in 
the IGM. 

We can also use the luminosity functions to estimate 
the fraction of UV continuum sources that have strong 
Lya emission. This fraction rises rapidly at low redshifts 
(from 5% at z = to ~ 20% at z = 2; Cowie, Barger & 
Hu 2010) but above z = 2 it has been less clear whether 
it continues to increase or drops again (Stark et al. 2010, 
2011; Zucca et al. 2009; Hayes et al. 2010). Above a 
UV continuum magnitude of —20.2 the Lya sample com- 
prises 24 ± 6% of the color-selected continuum sample 
at z = 6.5. For z = 5.7 the corresponding number is 
18 ± 2% above —19.9, where the errors are formal sta- 
tistical errors. These fractions are very similar to those 
found at z — 3 by Shapley et al. (2003) for sources with 
similar rest-frame Lya equivalent widths (above 20A), 
which would suggest very little high-redshift evolution. 
It would represent a drop from the values at z = 5.6 - 
6.1 found by Stark et al. (2011) who argue that nearly all 
continuum-selected objects at these redshifts are LAEs 
with rest-frame EW > 20A. However, we caution that 
the present numbers have substantial systematic uncer- 
tainties. The uncertainty in the UV to Lya conversion 
ratio can result in significant changes: thus if the true 
conversion was 0.03 (0.13), the second lowest (second 
highest) value observed in the sample, the normalization 
would rise (fall) to approximately 0.6 (0.1). In addition, 
as we have discussed above, other measurements of the 
Lya luminosity function at these redshifts (Kashikawa 
et al. 2006, 2011; Ouchi et al. 2008, 2010) have been a 
factor of two higher than Hu et al.'s and using these lu- 
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minosity functions would increase the LAE fraction by 
this amount. It is also possible that the UV continuum 
samples we are comparing with may still contain spurious 
objects, including red stars and lower redshift red galax- 
ies, which would also increase the LAE fraction, and that 
there may be substantial uncertainties owing to cosmic 
variance in these continuum luminosity functions. 

The present results show that LAEs at z = 6.5 have 
remarkably similar properties to those at lower redshifts. 
They are small, have moderately high star formation 
rates, and have substantially complete escape of the Lya 
photons. Other properties such as the shape and widths 
of the Lya lines also appear invariant (e.g. Ouchi et al. 
2010, Hu et al. 2010). The results also show that the 
continuum luminosities are comparable to the brightest 
values seen in color-selected continuum samples, suggest- 
ing that even the strongest star-forming galaxies at these 



redshifts may show Lya in emission. We also find that 
the LAEs are a significant fraction of the UV-continuum 
selected objects at these redshifts, similar, within the sys- 
tematic uncertainties, to values seen at lower redshifts. 
The observations show no sign that any effects associated 
with rcionization are being seen at z = 6.5. 

We gratefully acknowledge support from NSF grants 
AST-0709356 (L. L. C), AST-0687850 (E. M. LL), and 
AST-0607871 (A. S.) and from Hubble Space Telescope 
HST-GO-11108.01-A (E. M. H.) 



REFERENCES 



Bond, N. A., Gawiser, E., Gronwall, C, Ciardullo, R., Altmann, 

M., & Schawinski, K. 2009, ApJ, 705, 639. 
Bond, N. A., Feldmeier, J. J., Matkovic, A., Gronwall, C, 

Ciardullo, R., & Gawiser. E. 2010. ApJ . 716, L200. 
Bouwens, R., et al. 2010, arXiv: 1006.4360 
Cowie, L. L., Barger, A. J. & Hu, E. M. 2010, ApJ, 711, 928. 
Dcharveng, J.-M., et al. 2008, ApJ, 680, 1074. 
Chariot, S., & Fall, S. M. 1993, ApJ, 415, 580. 
Dayal, P., Maselli, A., &c Ferrara, A. 2011, MNRAS, 410, 830. 
Dijkstra, M. & Westra, E. 2010, MNRAS, 401, 2343. 
Dijkstra, M., & Wyithe, J. S. B. 2010, MNRAS, 408, 352 
Finkelstein, S. L. et al. 2007, ApJ, 660, 1023. 
Giavalisco, M., et al. 2004, ApJ, 600, L93. 

Hayes, M., Schaerer, D., Ostlin, G., Mas-Hesse, J. M., Atex, H., & 

Kunth, D. 2010, arXiv: 1010.4796 
Hu, E. M., Cowie, L. L, McMahon, R. G., Capak. P., Iwamuro, F., 

Kneib, J.-P., Maihara, T., & Motohara, K., 2002, ApJ, 568, L75 

(erratum 576, L99). 
Hu, E. M., Cowie, L. L., Barger, A. J., Capak, P., Kakazu, Y., & 

Trouille, L. 2010, ApJ, 725, 394. 
Kashikawa, N., et al. 2006, ApJ, 648, 7. 
Kashikawa, N., et al. 2011, arXiv:1104.2330 
Kennicutt, R. C, Jr. 1998, ARAA, 36, 189. 
Kunth, D. et al. 2003, ApJ, 597, 263. 

Laurscn, P., Sommcr-Larsen, J., & Razoumov, A. O. 2011, ApJ, 
728, 52. 



Lehnert, M. D. et al. 2010, Nature, 467, 940. 

Loeb, A., & Rybicki, G. B. 1999, AJ, 524, 527. 

Mackenty, J. W., Kimble, R. A., O'Connell, R. W., & Townsend, 

J. A. 2010, SPIE, 7731, 27. 
Neufeld, D. A. 1991, ApJ, 370, L85. 
Ostlin, G. et al. 2009, AJ, 138, 923. 
Ouchi, M, et al. 2008, ApJS, 176, 301. 
Ouchi, M, et al. 2009, ApJ, 706, 1136. 
Ouchi, M., et al. 2010, ApJ, in press. 
Schaerer, D., & Verhamme, A. 2008, A&A, 480, 369. 
Shapley, A. E., Steidel, C. C, Pettini, M., & Adelberger, K. L. 

2003, ApJ, 588, 65. 
Shimisaku, K. et al. 2006, PASJ, 58, 313. 
Stark, D. et al. 2010a, MNRAS, 408, 1628. 

Stark, D. P., Ellis, R. S. &; Ouchi, M. 2011, ApJ, 2011, 728, L2. 
Taniguchi, Y. et al. 2009, ApJ, 701, 915. 
Tilvi, V. et al. 2010, Ap J, 721, 1853. 
Vanzella, E. et al. 2010, larXiv: 1011. "55001 

Verhamme, A., Schaerer, D., & Maselli, A. 2006, A&A, 460, 397. 
Zheng, Z., Cen, R., Trac, H., & Miralda-Escude, J. 2010, ApJ, 716, 
574. 

Zucca, E., et al. 2009, A&A, 508, 1217. 



